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Abstract 

We investigate how models for neutrino-nucleus cross sections based on different 
assumptions for the nuclear dynamics affect the forecasted sensitivities to neutrino 
oscillation parameters at future neutrino facilities. We limit ourselves to the quasi- 
elastic regime, where the neutrino cross sections can be evaluated with less uncer- 
tainties, and discuss the sensitivity reach to #13 and 5 at a prototype low-7 /3-beam, 
mostly sensitive to the quasi-elastic regime. 
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1 Introduction 



In recent years the experimental study of neutrino oscillations has much contributed to 
our knowledge of particle physics by establishing non vanishing neutrino masses and by 
measuring or constraining the corresponding mixing angles. Within the domain of neutrino 
oscillations, the main goal of the next generation of facilities is the measurement of the 
mixing angle #13, which at present is only limited by an upper bound, and the observation 
of leptonic CP violation, for which we have no hints at the moment. These measurements 
are extremely challenging, since the smallness of the present bound on sin# 13 < 0.22 
at 3cr pQ constrains the signals to be very subleading. It is therefore necessary to have 
a good knowledge of the detector response function as well as an accurate knowledge 
of the beam to control the various systematic errors with a very good precision. For 
this task, novel v beams with extremely low backgrounds and systematic uncertainties 
have been proposed for future measurements based on the decay of muons (Neutrino 
factories [21 [3]) or (3 unstable ions (/3-Beams jl]). The v e fluxes from such decays would be 
extremely pure and the systematic errors very small, especially if the flux is normalized 
through a measurement with a near detector. Conversely, the systematic uncertainties 
associated with the detection process should be kept at the same accuracy level. It is 
usually argued that the f M /z/ e neutrino-nucleus cross section ratio can be understood at the 
level of some percent and that the v e cross section can be measured with similar precision 
at a near detector so that the signal can have systematic uncertainties of the order of 
a few percent at the far detector. This, however, does not take into account our present 
ignorance of the actual value of the cross section at low energies and how different models 
lead to forecasted sensitivities to the unknown parameters that differ by much larger 
margins than the few percent uncertainty normally considered. This question can severely 
impact the comparison of the relative performance of different facilities depending on the 
model adopted to parametrize their cross sections, particularly if the different facilities are 
sensitive to distinct energy regions. In this letter we want to discuss more in detail this 
question showing that, when a realistic model of nuclear dynamics is adopted, the neutrino 
cross section can sizably affect the forecasted precision measurement of #13 and the CP 
violating phase 5. It is beyond the scope of this letter to discuss (and critically revise) all 
possible models for neutrino interactions, so we restrict ourselves to the Relativistic Fermi 
Gas Model (RFG), variants of which are widely used in many MonteCarlo codes and, as 
examples of more refined calculations, to the one based on the Spectral Function approach 
(SF), on the Relativistic Mean Field (RMF) approximation and on the Random Phase 
Approximation (RPA), also including the effects of multinucleon contributions. The main 
features of these models are briefly summarized in Sect J2] whereas in Sect J3] we show how 
the measurement of #13 and 5 can be affected by the different descriptions of the nuclear 
dynamics. We draw our conclusions in SectJU 



2 Summary of the charged current neutrino-nucleus 
cross sections 

We work in the quasi-elastic regime (QE), which is of interest in many current and planned 
experiments (among them, MiniBooNE has already released its first cross section mea- 
surement |5j). The doubly-differential cross section, in which a neutrino carrying initial 
four- momentum k = (E u , k) scatters off a nuclear target to a state of four-momentum 



k' = (Eg, k ) can be written in Born approximation as follows: 



dQdEt 16 7T 2 |k| A 

where Gf is the Fermi constant and V u d is the CKM matrix element coupling u and d 
quarks. The leptonic tensor, that can be written in the form 

k^ k v + k v k^ 9fj,u(k • k ) i s^a^ k ^ k (2) 



is completely determined by lepton kinematics, whereas the nuclear tensor , containing 
all the information on strong interactions dynamics, describes the response of the target 
nucleus. Its definition involves the initial and final hadronic states |0) and \X), carrying 
four momenta po and px, respectively, as well as the nuclear electroweak current operator 



r A -. 



war = T,(^A\x)(x\r A \o)5^(po + q -p x ) , (3) 



X 



where the sum includes all hadronic final states. The calculation of W A V is a complicated 
task which deserves some approximation; quite often the impulse approximation (IA) 
scheme is adopted, based on the assumptions that at large enough q the target nucleus is 
seen by the probe as a collection of individual nucleons and that the particles produced 
at the interaction vertex and the recoiling (A — l)-nucleon system evolve independently. 
Within this picture, the nuclear current can be written as a sum of one-body currents, 
i.e. J A — > J2i Ji \ while the final state reduces to the direct product of the hadronic state 
produced at the weak vertex (with momentum p ) and that describing the (A — l)-nucleon 
residual system, with momentum p^: \X) — > \i,p ) <8> \71,Pk)- 



2.1 The Spectral Function approach 

The calculation of the weak tensor as described in Ref. [6] naturally leads to the concept 
of Spectral Function. In fact, the final expression of the hadronic tensor can be cast in 
the following form: 

WT = \ld'pdEP{p,E)—\ W^(p,q), (4) 



where E p = y |p| 2 + m 2 N and the function P(p, E) is the target Spectral Function, i.e. the 
probability distribution of finding a nucleon with momentum p and removal energy E in 
the target nucleus. It then encodes all the informations about the initial (struck) particle. 
The quantity W^ u is the tensor describing the weak interactions of the i-th nucleon in 
free space; the effect of nuclear binding of the struck nucleon is accounted for by the 
replacement q — (u, q) — y q — (p, q) with v = E\ p+ q\ — E\ p \. It follows that the second 
argument in the hadronic tensor is p = (E\ p \, p). Substituting Eq. @ into Eq. (JTJ, we get 
the final formula for the nuclear cross section: 

'tpiEPtoE)^, (5) 



dQdE e J r ' dQdE, 

in which we have redefined the elementary cross section as 

G^Oelem G\ V^ d |k | 1 



<KldE t 32 7r 2 |k| AE p E\ p+(i \ 



L^W^ . (6) 



The calculation of P(p,E) has been only carried out for A < 4 |7J; however, thanks to 
the simplifications associated with translation invariance, highly accurate results are also 
available for uniform nuclear matter, i.e. in the limit A — > oo with Z=A/2 [8] (Z denotes 
the number of protons). The Spectral Functions for medium- heavy nuclei have been 
modeled using the Local Density Approximation (LDA) [9J, in which the experimental 
information obtained from nucleon knock-out measurements is combined with the results 
of theoretical calculations of the nuclear matter P(p, E) at different densities. 

2.2 The Relativistic Fermi Gas 

The RFG JTU] model, widely used in MonteCarlo simulations, provides the simplest form 
of the Spectral Function: 



where pp is the Fermi momentum and Eb is the average binding energy, introduced to 
account for nuclear binding. The term in parenthesis is a constant needed to normalize 
the Spectral Function to the number of target nucleons, A. Thus, in this model pp and 
Eb are two parameters that are adjusted to reproduce the experimental data. For oxygen, 
the analysis of electron scattering data yields pp = 225 MeV and E B = 25 MeV [TT] . 

2.3 The Relativistic Mean Field approach 

Within the Relativistic Mean Field approximation we refer to the model described in [13] , 
where, like in the previous cases, the nuclear current is written as a sum of single-nucleon 
currents. The wave functions for the target and the residual nuclei are described in terms 
of an independent-particle model. Then, the transition matrix elements can be cast in the 
following form: 



where ip B and ipp are the wave functions for initial bound and final outgoing nucleons, 
respectively, and is the relativistic current operator. In particular, the relativistic 
bound-state wave functions (for both initial and outgoing nucleons) are obtained as a 
solution of the Dirac equation, in the presence of the same relativistic nuclear mean field 
potential, derived from a Lagrangian containing a, u and p mesons. The calculated cross 
sections correctly account for the inclusive cross section which are interested in. 

2.4 The Random Phase Approximation 

The last model we want to take into account has been introduced in [T5], where the 
hadronic tensor is expressed in terms of the nuclear response functions treated in the Ran- 
dom Phase Approximation (RPA). The response functions are related to the imaginary 
part of the corresponding full polarization propagators and the introduction of the RPA 
approximation means that the polarization propagators are the solutions of integral equa- 
tions involving the bare propagators and the effective interaction between particle-hole 
excitations. Within this formalism, the authors of [T5] were able to show that multinu- 
cleon terms sizably increase the genuine charged current QE cross section [151 02] m such 







a way to reproduce the MiniBooNE results [5]. The mechanism responsible for the en- 
hancement that brings the theoretical cross section into agreement with the data is multi 
nucleon knock out, leading to two particle-two hole (2p2h) nuclear final states. In the 
following, we will refer to this "generalized" QE cross section as RPA-2p2h whereas we 
adopt the short RPA for the genuine QE cross section. 

2.5 Comparison of the cross sections 

To summarize this section, we present in FigJT] a comparison of the five total QE cross 
sections for the 16 — > fi~X process (left panel) and 16 — > fi + X (right panel), in 
the energy range E v ~ [0,0.75] GeV. The curves have been computed using the dipole 
structure of the form factors and, in particular, a value of the axial mass close to rriA ~ 1 
GeV. As it can be easily seen in the left panel, the RFG prediction sizably overestimates 
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Figure 1: Total charged current QE cross sections for the ujfO — > fi~X (left panel) and 
ujfO — » fi + X (right panel) processes in the energy range E v ~ [0, 0.75] GeV. 

the SF, RMF and RPA results by roughly 15%, a fact that is well known to happen also for 
many other models with a more accurate description of the nuclear dynamics than the RFG 
approach (see, for instance, [14]). On the other hand, the inclusion of 2p2h contributions 
largely enhances the QE cross section in the RPA approximation for energies above ~ 0.5 
GeV, although it is smaller than the RFG at smaller energies. For antineutrinos (right 
panel) the observed pattern is almost the same. From this comparison, a qualitative 
understanding of the impact of different models of neutrino-nucleus cross sections on the 
forecasted sensitivity of future facilities can be already derived. A quantitative analysis of 
these effects will be the subject of the next section. 

3 The impact on the (9is-S) measurement 

To estimate the impact of different models of the cross section on the measurement of 
#i3 and leptonic CP violation we choose a 7 = 100 /3-Beam facility as a representative 
example. The choice is motivated because the neutrino flux from such a facility spans up 
to ~ 0.7 GeV with the peak around 0.3 GeV and is thus mostly sensitive to the quasielastic 
region explored here. The /3-Beam concept was first introduced in Ref. [4] and involves the 
production of /3-unstable ions, accelerating them to some reference energy, and allowing 



them to decay in the straight section of a storage ring, resulting in a very intense and 
pure u e or v e beams. The v e — > "golden channel", which has been identified as the 



most sensitive to all the unknown parameters [18], can be probed at a far detector. We 
have considered here the original /3-Beam proposal, where 18 Ne ( 6 He) ions are accelerated 
to 7 ~ 100 at the CERN SPS and stored so that v e (z/ e ) beams are produced and the 
golden channel oscillation is searched for at a Mton class water Cerenkov detector located 
at L — 130 km at the Frejus site, detailed analyses of the physics performance of this 
setup can be found in Refs. [19l [2U1 [2H E21 [231 12H [2S1 [261 123 [28]. In order to simulate 
the detector response when exposed to such a beam both in terms of signal efficiency and 
background, we have used the migration matrices derived in Ref. [29]. Systematic errors 
of 2.5% and 5% in the signal and background respectively have been taken into account. 
In all the simulations the following best fit values and la errors for the known oscillation 
parameters were assumed pQ Am^ = (7.6 ± 0.2) • 10~ 5 eV 2 , Am| x = (2.5 ± 0.1) • 10~ 3 
eV 2 , 9 12 = 34.0 ± 1.0 and #23 = 45.0 ± 3.6. These parameters were marginalized over to 
present the final curves. The evaluation of the performance of the facility made use of the 
GLoBES software [3"0"II3"T| . It is important to notice that we are only using the quasielastic 
contribution to the neutrino cross section depicted in Fig{TJ 

As an illustration we have focused on the dependence on the nuclear model adopted 
of two different observables, namely the CP and 9i 3 discovery potentials, defined as the 
values of the CP- violating phase Sep and #13 for which respectively the hypothesis of CP 
conservation Sep = 0, ±7r or # 13 = can be excluded at 3cr after marginalizing over all 
other parameters. The CP discovery potential is shown in FigfJJ where we superimposed 
the results obtained using the RFG cross section and the SF, RMF and RPA calculations. 
In the left panel, the CP discovery potential is represented in the (6*13, £c*p)-plane; we 
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Figure 2: Left panel: CP discovery potential in the {Q\3,5cp) plane. Solid lines refer to 
the SF model, dotted lines to the RFG, short-dashed lines to the RMF, dot-dashed to RPA 
and long-dashed to RPA-2p2h. Right panel: the CP fraction. 

clearly see that, for 5c p ~ ±90° (where the sensitivity is maximal) the RFG model gives a 
prediction which is around a factor 2 better than the SF, RMF and RPA models (which, 
as expected, behave almost in the same way) in sin 2 26*13 an d around a 40% better than 
the RPA-2p2h. This is not surprising because the /3-Beam facility used in our simulations 
mainly probes energies smaller than 0.5 GeV, where the RFG is still larger than any other 



model (see Fig{T]). For the other points in the parameter space, the difference is less evident 
but still significant. The same information can be summarized in the right panel making 
use of the 5 fraction (Sf), that represents the fraction of values of Sep for which CP can 
be discovered at a given # 13 . The 5 fraction has a maximum around sin 2 2# 13 ~ 10 -2 
where the RFG model predicts 5p ~ 0.75 and the other models give 5p ~ 0.7. The 
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Figure 3: Left panel: 9i 3 discovery potential in the (9i 3 ,5cp)-plane. Solid lines refer to 
the SF model, dotted lines to the RFG, short-dashed lines to the RMF, dot-dashed to RPA 
and long-dashed to RPA-2p2h. Right panel: the CP fraction. 

difference in the sensitivity to #13 can be seen in Fig. [3] where we show the results in 
the (sin 2 2# 13 , 5<7p)-plane. In this case the predictions of the SF, RMF and RPA models 
differ by up to a factor of ~ 60% compared to the RFG model for 5c p ~ ±90° while the 
difference is less pronounced for 5c p ~ 0°. The RPA-2p2h results are more similar to the 
RFG for 5cp ~ 0° and to the other models for 5cp ~ 180°. This is also evident in the 
right panel where we show the CP fraction. 

Finally, it is interesting to observe the effect of using different nuclear models also in 
the simultaneous determination of #13 and 5c p- We present an example in FigJH where, 
for the input value (#13, 5cp) = (0.9°, 30°), indicated with a dot, we show the capability of 
the /3-beam to reconstruct the true values of our observables at 90% CL. For the sake of 
simplicity, we did not include degenerate solutions coming from our ignorance of the octant 
of #2 3 and the hierarchy in the neutrino mass ordering but only the so-called "intrinsic" 
one. The main feature here is that using the RFG and RPA-2p2h models we are able to 
reconstruct the true values of #13 and 5c p within reasonable uncertainties, whereas with the 
other models we can only measure two distinct disconnected regions (the fake one around 
the value of #13 and 5cp ~ 180°), which worsen the global sensitivity on those parameters. 
The effects we have mentioned have been generalized to account for the following cases: 

• the value of the axial mass does not fill the gap among the SF approach and the 
RFG model, at least in the range e [1, 1.2] GeV. This is not surprising and has 
been already the subject of an extensive analysis [P7] ; 

• the same sensitivity behaviour as observed in FigsfSJIHlhas been also seen for a differ- 
ent nuclear target, namely 56 Fe. This points to the conclusion that the RFG model 
overestimates the sensitivities to 6*13 and 5cp for a vast class of nuclear targets. 
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Figure 4: 90% CL contour for the input value (613, 5cp) = (0.9°, 30°). Solid lines refer to 
the SF model, dotted lines to the RFG, short-dashed lines to the RMF, dot-dashed to RPA 
and long-dashed to RPA-2p2h. 

4 Conclusion 

In this letter we have analyzed the impact of different neutrino-nucleus charged current 
QE cross sections on the forecasted sensitivity of the future neutrino facilities to the 
parameters # 13 and b~cp- We considered five different calculations, based on the RFG model 
(widely used in the MonteCarlo codes) and those based on the Spectral Function approach, 
on the Relativistic Mean Field approximation and of the Random Phase Approximation 
(including multinucleon contributions, also) in the quasi-elastic regime. We found that the 
sensitivities computed from SF, RMF and RPA (without the multinucleon contribution) 
models are worse than the RFG results, for both 13 and Sep, by up to a factor 2 when 
an Oxygen nuclear target is used to compute the event rates. To a less extent, this is 
also true for the RPA-2p2h model, due to the fact that the quasi-elastic cross section 
is larger than the genuine quasi-elastic cross section. These variations with the nuclear 
models are so large that they cannot be taken into account with the few percent systematic 
error uncertainty expected from dedicated measurements at near detectors. Indeed, the 
present uncertainty is much larger and allows for discrepancies between models that can 
strongly affect the forecasted sensitivity of a given facility, implying that special care must 
be adopted when comparing the performance of different facilities. We have also checked 
that other nuclear targets produce quantitatively similar results. This suggests that the 
use in MonteCarlo simulations of more refined nuclear models for the neutrino-nucleus 
interaction is mandatory, especially for those facilities whose bulk of events is in the quasi- 
elastic regime. 
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